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Abstract 

This paper discusses a FORTRAN program that 
provides the user with a powerful and flexible 
tool for the linearization of aircraft models. 

The program LINEAR numerically determines a linear 
systems model using nonlinear equations of motion 
and a user-supplied, nonlinear aerodynamic model. 
The system model determined by LINEAR consists of 
matrices for both the state and observation equa- 
tions. The program has been designed tu allow 
easy selection and definition of the state, con- 
trol, and observation variables to be used In a 
particular model. Also, Included in the report is 
a comparison of linear and nonlinear models for a 
high-performance aircraft. 

Symbol s 

The units associated with the quantities 
listed below are expressed in a generalized 
system. Where applicable, quantities are defined 
with respect to the body-, xis system. 

Variables 


A state matrix of the equation, 

x * Ax ♦ Bu 

A axial force, lb 

A' state matrix of the equation, 

Cx * A'x ♦ B'u 

a n> i normal acceleration not at vehicle 
center of gravity, g 

a n x-body axis accelerometer at vehicle 

center of gravity, g 

a„ x • x-body axis accelerometer not at vehicle 
center of gravity, g 

a n y-body axis accelerometer at vehicle 

y center of gravity, g 

a n . y-body axis accelerometer not at 
y * vehicle center of g r av1ty, g 

a nz z-body axis accelerometer at vehicle 

center of gravity, g 

an z ( z-body axis accelerometer not at 
vehicle center of gravity, g 

•Aerospace Engineer. Member AIAA. 

tAerospace Engineer. 


a x acceleration along the x-body axis, g 

a y acceleration along the y-body axis, g 

a 2 acceleration along the z-body axis, g 

B control matrix of the equation, 

x * Ax + Bu 

B' control matrix of the equation, 

Cx = A'x + B'u 

C C-matrlx of the state equation 

Cx * A'x + B'u 

D dynamic interaction matrix for state 

equation, x » Ax + Bu ♦ Dv 

D' dynamic interaction matrix for the 

state equation, Cx = A'x + B'u + D'v 

D drag force, lb 

E dynamic interaction matrix for the 

observation equation, 
y = Hx ♦ Fu ♦ Ev 

E' dynamic interaction matrix for the 

observation equation, 

y » H'x + Gx + F'u + E'v 

E s specific energy, ft 

F feed-forward matrix of the equation, 

y = Hx + Fu 

F' feed-forward matrix of the observation 
equation, y * H'x + Gx* + F'u 

f nonlinear state vector function 

fpa flightpath acceleration, g 

G G-matrlx of the observation equation, 

y ■ H'x + Gx ♦ F'u 

g nonlinear observation vector function 

H observation matrix of the equation, 

y = Hx ♦ Fu 

H' observation matrix of the equation, 

y » H'x + Gx’ + F'u 
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h altitude, ft 

h,1 altitude not at vehicle center of 

gravity, ft 

L lift, lb 

M Mach number 

N normal force, lb 

n load factor 

P $ specific power, ft/sec 

p roll rate, deg 

P a ambient pressure, 1 b/f t 2 

P s stability axis roll rate, deg/sec 

P t total pressure, lb/ft? 

q pitch rate, deg/sec 

q dynamic pressure, lb/ft? 

q c impact pressure, 1 0 / f 1 2 

q 5 stability axis pitch rate, deg/sec 

Re Reynolds number 

Re" Reynolds number per unit length, ft'l 

r yaw rate, deg/sec 

r s stability axis yaw rate, deg/sec 

T ambient temperature, °R 

T total angular momentum, slug-ft2/sf.x2 

T t total temperature, °R 

t time, sec 

u velocity In x-axIs direction, ft/sec 

u control vector 

u 0 control vector at analysis point 

V total velocity, ft/sec 

v dynamic interaction vector 

v velocity In y-axls direction, ft/sec 

w velocity in z-axls direction, ft/sec 

X total force along the x body axis, lb 

x state vector 


x 0 state vector along nominal trajectory 

Y side force, lb 

y observation vector 

y 0 observation vector along nominal trajectory 

Z total force along the z-body axis, lb 

a angle of attack, deg 

а, 1 angle of attack measurement not at 

vehicle center of gravity, deg 

0 angle of sideslip, radians 

б, 1 angle of sideslip measurement not at 

vehicle center of gravity, deg 

Y flight path angle, dog 

A small perturbation value 

«L Incremental rolling moment, ft lb 
6M Incremental pitching moment, ft lb 

6N Incremental yawing moment, ft lb 

6X Incremental x-body axis force, lb 

6Y Incremental y-body axis force, lb 

61 Incremental z-body axis force, lb 

e pitch angle, deg 

+ roll angle, deg 

heading angle, deg 


Superscripts 

. Indicates derivative with respect to time 

T transpose of a vector or matrix 


Introduction 

The program LINEAR, described In this paper, 
was developed at the Dryden Flight Research 
Facility o f the NASA Ames Research Center to pro- 
vide a standard, documented, and verified tool to 
derive linear models for aircraft stability analy- 
sis and control 1 ’w design. This development was 
undertaken to eliminate the need for the aircraft- 
specific linearization programs common In the aero- 
space Industry. The lack of available, documented 
linearization program provided a strong motivation 
for this paper and a NASA technical paper (to be 
published) on the deve’opment of LINEAR. Refer- 
ence 1 represents the only available documented 
linearization program before LINEAR. 
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Linear system models of aircraft dynamics and 
sensors are an essential part of both vehicle- 
stability analysis and control law design. These 
models define the aircraft system in the neighbor- 
hood of an analysis point and are determined by 
the linearization of the nonlinear equations 
defining vehicle dynamics and sensors. This 
report describes a FORTRAN program that provides 
the user with a powerful and flexible tool for the 
linearization of aircraft models. LINEAR is a 
program with well-defined and generalized inter- 
faces to aerodynamic and engine models and is 
designed to address a wide range of problems with- 
out the requirement of program modification. 

The program LINEAR numerically determines a 
linear systems model using nonlinear equations of 
motion and a user-suppl iej nonlinear aerodynamic 
model. lINEAR Is also capable of extracting both 
linearized engine effects (such as net thrust, 
torque, and gyroscopic effects) and including 
these effects in the linear system model. The 
point at which this linear system model is defined 
is determined either by completely specifying the 
state and control variables or by specifying an 
analysis point on a trajectory, selecting a trim 
option, and directing the program to determine the 
control variables and remaining state variables. 

The system model determined by LINEAR consists 
of matrices for both the state and observation 
equations. The program has been designed to pro- 
vide an easy selection and definition of the 
state, control, and observation variables to be 
used in a particular model. Thus, the order of 
the system model is completely under the user's 
control. Further, the program provides the flexi- 
bility of allowing alternate formulations of both 
the state and observation equations. 


Program Overview 

The program LINEAR numerically oetermines a 
linear systems model using nonlinear equations of 
motion and a user-supplied nonlinear aerodynamic 
model (Fig. 1). LINEAR is also capable of extract- 
ing linearized gross-engine effects, (such as net 
thrust, torque, and gyroscopic effects) and Includ- 
ing these effects in the linear system model. The 
point at which this linear system model is defined 
is determined either by specifying the state and 
control variables or by selecting an analysis point 
on a trajectory, selecting a trim option, and 
allowing the program to determine the control vari- 
ables and remaining state variables to satisify the 
trim option selected. 

Because the program is designed to statisfy the 
needs of a broad class of users, a wide variety of 
options has been provided. Perhaps the most impor- 
tant of these options are those that allow user spe- 
cification of the state, control, and observation 
variables to be included in the model derived by 
LINEAR. 

Within the program, the nonlinear equations of 
motion include 12 states representing a rigid air- 
craft flying in a stationary atmosphere over a flat, 
nonrotatiny earth. Internally, the state vector x 
is computed as 

* * [p. q. r, V, a, B, ♦, 8, i|>, h, x, y] T 

The internal -control vector, u, can contain up 
to 30 controls. The internal observation vector, y, 
contains 120 variables including the state vari- 
ables, the time derivatives of the state variables, 
the control variables, and a variety of other para- 
meters of interest. Within the program. 


LINEAR has several features that make it 
unique among the standard linearization programs 
in the aerospace industry. The most significant 
of these features is flexibility. By generalizing 
the surface definitions and making no assumptions 
of symmetric mass distributions, the program can 
be applied to any aircraft in any phase of flight, 
except hover. The unique trimming capababi 1 ity , 
provided by means of a user supplied subroutine, 
allow unlimited possibilities of trimming strate- 
gies and surface scheduling, which are particu- 
larly important for oblique-winged vehicles and 
aircraft having multiple surfaces effecting a 
single axis. The formulation of the equations of 
motion permit the inclusion of thrust vectoring 
effects. The ability to select, without program 
modification, the state, control, and observation 
variables for the linear models, which when com- 
bined with the large number of observation quan- 
tities available, allows any analysis problem to 
be attacked with ease. 

This paper provides an introduction to the 
program LINEAR. The trimming capabilities of 
LINEAR are discussed frem both a theoretical and 
implementation perspective. Time history com- 
parisons of linear and nonlinear models of the 
high-performance aircraft were used as part of 
the validation of LINEAR and are presented to 
illustrate the capabilities o : the program. 



where 


x • a y» a n x * a n y . a n z « a n> a n x j» a n y> i» 
a n z ,i . a n,i • nJ T 

y 2 = [a. Re, R£, M, q, q c , p a , q c /p a , p t , T, T t , 

Ve. V C ] T 


y 3 = [y. fpa, y , h, h/57.3] 
y 4 = [E s . P S ] T 
y 5 * [L. 0 , N, A] T 
y 6 * [u, v, w, u, v, w] T 

• J 

y 7 = [<X,1, B, i, h,i, h,i] 
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and 


y 6 * [T. p s . q s .r s ] T 

From the Internal formulation of the state, 
control, and observation variables, the user must 
select the specific variables desired in the out- 
put linear model. Figure 2 illustrates the selec- 
tion of the variables in the state vector for a 
requested linear model. From the internal formu- 
lation on the right, the requested model is con- 
structed, and the linear system matrices are 
selected in accordance with the user specification 
of the state, control, and observation variables. 

The model derived by LINEAR is determined at 
an analysis ;,oint. LINEAR allows this analysis 
point to be defined as a true, steady-state condi- 
tion on a specified trajectory (a point at which 
the rotational and translational accelerations are 
zero) or a totally arbitrary state on an arbitrary 
trajectory. The program LINEAR provides the user 
with several options described in detail in the 
"Analysis-Point Definition" section of this paper. 
These analysis-point-definition options allow the 
user to trim the aircraft in "wings-1 evel " flight, 
pushovers, pullups, level turns, or zero-sideslip 
maneuvers; also included is a nonfrimming option in 
which the user defines a totally arbitrary condi- 
tion about which the linear model is to be derived. 

The linear system matrices are determined bv 
numerical perturbation and are the first-order 
terms of a Taylor series expansion about the analy- 
sis point as described in the "Linear Models" sec- 
tion of this paper. The formulation of the output- 
system model is under user control. The user can 
select state-equation matrices corresponding to 
either the standard formulation of the state 
equation, 

x = Ax ♦ Bu 

or the generalized equation, 

Cx = A'x + B'u 

The observation matrices can be selected from 
either of two formulations corresponding to the 
standard equation, 

y = Hx + Fu 

or the generalized equation, 

y * H'x ♦ Gx + F'u 

In addition to the linear system matrices, 
LINEAR also computes the nondimensional stability 
and control derivatives at the analysis point. 

The input file for the batch version of LINEAR 
is a Curd-image file tnat defines the geometry 
and mass properties of the aircraft and selects 
various program options. Within this input file 
the state, control, and observation vectors desired 
in the output linear model are defined and the 


analysis-point options are selected. The interac- 
tive version of LINEAR allows user Inputs from a 
terminal, the use of a card-image input file, or a 
combination of both. 

The output of LINEAR is two files: one con- 

taining the linear system matrices and one docu- 
menting the options and analysis points selected 
by the user. The former of these files is 
intended to be used with follow-on design and 
analysis programs. The latter of these files con- 
tains all the information provided on the former 
file and also includes the details of the analysis 
point, the nondimensional stability, and control 
derivatives. 

To execute LINEAR, four user-supplied subrou- 
tines are required. These subroutines define the 
nonlinear aerodynamic model, the gross-engine 
model, and the gearing between the LINEAR trim 
inputs and the surfaces modeled in the aerodynamic 
model. The gearing model. Illustrated in Fig. 3, 
defines how the LINEAR trim inputs will be con- 
nected to the surface models and allows schedules 
and nonstandard trimming schemes to be employed. 

Li near Models 

The linearized system matrices computed by 
LINEAR are the first-order terms of a Taylor 
series expansion about the analysis point^.3 and 
are assumed to result in a time-invariant linear 
system. The validity of this assumption is dis- 
cussed in the "Ana’ysis-Point Definition" section 
of this paper. The technique employed to obtain 
these matrices numerically is a simple approxima- 
tion to the partial derivative, 

3f _ f(x 0 + Ax) - f(x 0 - Ax) 

3 X " Fax' ' ~ 

where f is a general function of x, an arbitrary, 
independent variable. 

From the generalized nonlinear state, 

Tx * f (x, x, u) 
and observation “quations, 

y * g (x, x, u) 

the program determines the linearized matrices for 
the generalized formulation of the system: 

C6x = A'6x + B'6u 
Sy a H' 6x G 6x + F ‘ 6u 

where 




A' 


3f 
* 7x 


B' 



H' 


-ia 

3x 


with all derivatives evaluated along the nomi- 
nal trajectory defined by the analysis point 

( *o • *0 • u o ) • 

By determining the matrices of the generalized 
formulation of the system, LINEAR avoids the pro- 
blems Inherent In systems having the nonlinear 
formulation 


x * f ( x , x\ u) 


G 


. la 

3l 


and 


r 


3g 

3u 


with all derivatives evaluated along the nomi- 
nal trajectory defined by the analysis point 

(*o. *o, u 0 ) and the state, time derivative of 
state, and control vectors can be expressed as 
small perturbations about the nominal trajectory, 
so that. 


y * g(x , x, u) 

rather than 

x* = f(x, u) 
y = g(x, u) 

The latter formulation is often assumed implic- 
itly, although for most aircraft the former formu- 
lation Is correct particularly because of the pre- 
sence of the aerodynamic derivatives with respect 

to angle-of-attack rate (a) and angle of sideslip 

rate (6). If the system matrices of the standard 
formulation of the linear system, given by 


r. ~ x 0 *■ fix 


x = Ax + Bu 


• • • 
X * x 0 + fix 


and 


u * u 0 + fiu 


In addition to the matrices for the generalized 
system just described, the user has the option of 
requesting linearized matrices for a standard for- 
mulation of the system: 


fix = A fix + B 6u 
fiy = H fix + F fiu 


where 


and 


T - 


3f 


3x 


•1 3f 
ix 


B 


3f -1 
3x 


3f 

3u 




sin - 1 

3x 


3f 

3x 




_3f - 1 
3x 


3f 

3u 


y = Hx + Fu 

are determined directly from nonlinear functions 

having state rate (x) dependence, an iterative 
scheme must be employed with each perturbation of 
the state and control variables while computing 
the partial derivatives. This is not always done, 
making the resultant linear system matrices 
incorrect. 

LINEAR also provides two nonstandard matrices, 
D and E, in the equations, 

x = Ax + Bu + D v 
y = Hx + Fu + E v 
or, D'and E', in the equations, 

Cx = A'x + D'u + D'v 


y = H'x + Gx + F'u + E'v 

These dynamic interaction matrices include the 
effect of external forces and moments acting on 
the vehicle. The components of the dynamic 
interaction vector, v, are incremental body-axis 
forces and moments: 

rsi 


6M 

6N 
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Thus 


O'- 


3f 

37 


and 



D 


T 


3f -1 
3x 


3f 

■57 


and 


E 


ia ♦ ia 

3V 3x 


T 


3f 
3 x 



The purpose of these matrices is to allow the 
effects of unusual subsystems or control effectors 
to be easily Included In the vehicle dynamics. 


nonequilibrium condition raises theoretical Issues 
beyond the scope of this report, but of which the 
potential user should be awa r e. While all of the 
analysis-point definition options provided In 
LINEAR have been found to be useful In the analysis 
of vehicle dynamics, not all of the linear models 
derived about these analysis points result In the 
time-invariant systems assumed In this report. 
However, the results of the linearization provided 
by LINEAR do give the appearance of being time 
Invariant. 

The linearization prccess as defined In this 
report Is always valid for some time Interval 
beyond the point In the trajectory about which the 
linearization is done. However, for the resultant 
system to be truly time Invariant the vehicle must 
be in a sustainable, steady-state flight condition. 
This requirement Is something more than merely a 
trim requirement, which Is typically represented as 

x(t) - 0 


The default output matrices for LINEAR are 
those for the standard system formulation. 

However, the user can select matrices for either 
the generalized or the standard state and obser- 
vation equations In any combination. Internally, 
the matrices are computed for the generalized 
system formjlation and then combined appropriately 
to accommodate the sysiem formulation requested by 
the user. 

Analysis-Point Definition 

The point at which the nonlinear system 
equations are linearized is referred to as the 
analysis point. This point can represent a true, 
steady-state condition on the specified trajectory 
(a point at which the rotational and translational 

acceleration are zero) 4 * 5 or a totally arbitrary 
state on a trajectory. LINEAR allows the user to 
select from a variety of analysis points. Within 
the program, these analysis points are referred to 
as trim conditions, and several trim options are 
available to the user. The arbitrary state and 
control option Is designated NOTRIM, and in select- 
ing this option the user must specify all nonzero 
state and control variables. For the equilibrium 
conditions, the user specifies a minimum number of 
parameters and the program numerically determines 
required state and control variables to force the 
rotational and translational accelerations tj zero. 
The analysis-point options are described In detail 
in the following subsections. 


Indicating that for trim, all the time derivatives 
of the state achieved must be zero. (This Is not 
the case, however: Trim Is achieved when the 

acceleration-like terms are Identically zero; no 
constraints need to be placed on the velocity-like 

terms In x. Thus, for the model used in L I NEAP , 

only p, q, r, tf, a, and § must be zero to satisfy 
the trim condition.) The trim condition is 
achieved for the straight and level, pushover/ 
pullup, level turn, thrust-stabilized turn, and the 
beta-trim options listed below. In general, the 
analysis-point-definition options for NOTRIM and 
specific-power conditions do not result In a trim 
condition. 

Of these analysis-point options resulting in a 
trim condition, only the straight-and-level and 
level -turn options force the model to represent 
sustainable flight conditions. In fact, only In 
the special case where the fllghtpath angle Is zero 
does this occur for these options. 

As stated earlier, the linearization of a 
nonlinear model and its representation as a time- 
invariant system is always valid for some time 
interval beyond the analysis point on the trajec- 
tory. This time interval is determined by several 
factors such as trim, sustainable flight condi- 
tions, and in the end by accuracy requirements 
placed on the representation. Thus, in using the 
linear models provided by this program, the user 
should exercise some caution. 


For all of the analysis-point definition 
options, any state or control parameter may be 
Input by the user. Those state variables not 
required to define the analysis point are used as 
initial estimates for the calculation of the state 
and control conditions that result in zero rota- 
tional and translational accelerations. 

It should be noted that the option of allowing 
the user to linearize the system equations about a 


Untrimmed 

For the untrimmed option, the user specifies 
all state and control variables that are to be set 
at some value other than zero. The number of state 
variables specified Is entirely at the user's dis- 
cretion. If any of the control variables are to 
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be nonzero, the user must specify the control rate. The two options available allow the user to 

parameter and Its value. The untrlmmed option specify either the angle of attack or the load 

allows the user to analyze the vehicle dynamics factor for the analysis point. The altitude and 

at any flight condition, Including transitory Mach number at the analysis point must be spe- 

condltlons. clfled for either of the two options. The user 

also must specify the value of the thrust-trim 
Stralght- a nd-Level Trim parameter. 


The s :raight-and-level trims available In 
LINEAR are constant fl Ightpath-angle trims at 
"wlngs-level ." Both option;, available for 
stral ght-and-level trim allow the user to specify 
either a fllghtpath angle or an altitude rate. 
However, since the default value frr these t*rms 
Is zero the default for both types ?f stra ^ht- 
and-level trim Is wlngs-level, hoi iz ntal flight. 

The two options available for stralght-and- 
level trim require the user to specify altitude 
and either an angle of attoe* nr Mach number. If 
a specific angle-of-attack and altU^-'* combina- 
tion Is desired LINEAR determines the velocity 
required for the requested trajectory. LINEAR 
also takes the user's specified Mach number and 
altitude, and uses the trim routines In the 
program to determine the angle of attack needed 
for the requested trajectory. 

Pushover/Pu 1 1 up 

The analysis-point -definition options for 
pushover/pul lup result In wings-level flight at 
load factors of greater than or less than one. 

For the elevated-load-factor case, the analysis 
point is the minimum altitude point of a pullup 
when altitude rate Is zero. For the case of load 
factor less than one, this point results In a 
pushover at the maximum altitude with an altitude 
rate of zero. There are two options available for 
the analysis-point -definition of pushover/pul 1 up: 
(1) in which angle of attack is determined from 
the specified altitude and Mach number, and (2) in 
which angle of attack, altitude, and Mach number 
are specified and load factor is determined 
according to the constraint equations. 

Leve l Turn 

The analysis-point-definition options for 
level turn result in “nonwings-level constant- 
turn-rate flight at load factors greater than one. 
T.ie vehicle model is assumed to have sufficient 
excess thrust to trim at the condition specified. 
If thrust is not sufficient, trim will not result 
and the analysis point thus defined will have a 
nonzero (in fact, negative) velocity rate. 

The level trim options available ir LINEAR 
require the specification of an altitude and Mac.i 
number. The user can then use either angle of 
attack or load factor to define the desired flight 
condition. For either of these options the user 
may also request a specific fllghtpath angle or 
altitude rate. Thus, these analysis-point defini- 
tions may result in ascending or descending 
spirals, although the default is far the constant 
altitude turn. 


Beta Trim 


The anolysls-polnt-deflnltion option for beta 
trim results In "nonwings-level," horizontal flight 
with a zero-heading rate at a user-specified Mach 
number, altitude, and angle of sideslip. This 
trim option is nominally at 1 g, but as beta 
varies from zero, normal acceleration decreases 
and lateral acceleration Increases. For an aero- 
dynamically symmetric aircraft, a trim-to-zero 
beta using the heta-trlm option results In the 
same trimmed condition as the straight-and-level 
trim. However, for an asymmetric aircraft such as 
an oblique-winged vehicle, the t.*o trim options 
are not equivalent. 


Specific Power 


The analysis-point-definition option for spe- 
cific power results in a level turn at a user- 
specified Mach number, altitude, thrust-trim 
parameter, and specific power. Unlike tiie other 
trim options provided in LINEAR, the specific 
power option in general does not attempt to 
achieve a zero velocity rate. Because the alti- 
tude rate is zero and specific power is defined by 



the resultant velocity rate, V, will be 



However, the other acceleration-like ‘erms ( p , q, 

r, a, and B) will be zero if the requested analy- 
sis point is achieved. 


Program Validation 


Two methods were used to validate LINEAR. The 
first of these methods required that the elements 
of the matrices of the generalized, linear system 
equations be derived analytically. These results 
were compared to the numerically derived matrix 
elements produced by LINEAR. The second method of 
validation was by comparison of linear and nonli- 
near time histories. The linear time histories 
were generated by solving the linear system dif- 
ferential equations using 



and 


Thrust -vStabj li zed Turn 

The analysis-point -def in ition for a thrust- 
stabilized turn results in a constant -throttle 
"nonwings-level" turn with a nonzero altitude 


y(t) - y 0 + He At «x 0 ♦ f He A(t ' T) B 6u(t) di 
Jo 


♦ F 6u(t) 


jn 
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The nonlinear time histories were generated using 
Independent Implementations of the equations of 
motion and sensor models. However, the nonlinear 
simulations and the linearized matrices were based 
on the same Implementation of the nonlinear aerody- 
namic models. The nonlinear simulations used a 
12-state i.odel. The linear simulations used a 
nine-state model excluding T, x, and y. 

The simulation used for the time history com- 
parisons was based on a model of a two-engine, 
high-performance aircraft capable of speeds In 
excess of 2.0 Mach number and altitudes in excess 
of 50,000 ft. Two sets of time histories are pre- 
sented In Figs. 4 and 5. The first time history 
Is based on a stralght-and-level trim at 0.7 Mach 
number at 20,000-ft altitude. Figure 4 includes a 
time history from a four-state linear model In 
addition to the nonlinear and nine-state linear 
models. For the time history shown In Fig. 4 this 
model was 

m 


LeJ 

Figure 4 shows the results of a *2° elevator 
doublet. The first feature to note about the time 
histories In Fig. 4 Is that except for the velo- 
city and altitude time histories, there Is very 
little difference between the response of the 
three models. The second feature to note aboui 
these time histories Is that, except for the alti- 
tude response, the performance of the nlre-state 
and four-state models give Identical results. 

Figure 5 is based on a 3-g level turn at 0.9 Mach 
number at 35,000-ft altitude with a t 2° aileron 
doublet. Here, the time histories lor the lateral - 
directional parameter are almost identical, except 
for the bank angle; the longitudinal -parameter 
time histories, while in reasonable agreement, 
show steady-state divergence, but with similar 
transient response. 

The comparison of linear and nonlinear models 
in Figs. 4 and 5 illustrate two of the trim options 
available in LINEAR. These time histories also 
demonstrate the ability of LINEAR to derive linear 
aircraft models. The lack of complete agreement 
between the linear and nonlinear reflects on the 
adequacy of using linear models ratner than on 
LINEAR. If linear models are desired, the program 
LINEAR provides a useful tool for generating them. 

Concluding Remarks 

The FORTRAN program, LINEAR, was developed to 
provide a flexible, powerful, and documented tool 
to derive linear models for aircraft -stability 
analysis and control law design. The program 
LINEAR numerically determines a linear systems 
model using nonlinear equations of motion and a 
user-supplied, nonlinear, aerodynamic model. 

LINEAR is also capable of extracting both linear- 
ized engine effects (such as net thrust, torque. 


and gyroscopic effects) and Including these effects 
In tne linear system model. The point at which 
this linear system model Is defined Is determined 
either by completely specifying the state and con- 
trol variables or by specifying an analysis point 
on a tiajectory, selecting a trim option, and 
directing the program to determine the control 
variables and remaining state variables. 

The system model determined by LINEAR consists 
of natrlces for both the state and observation 
equations. The program has been designed to pro- 
vide an easy select'on and definition of the state, 
control, and observation variables to be used In a 
particular .'Odel. Thus, the order of the system 
model Is completely under user control. Further, 
the program provides the flexibility of allowing 
alternate formulations of both the state and obser- 
vation equations. 

LINEAR has several Matures that make It unique 
among the linearization programs common In the 
aerospace Industry. The most significant of these 
features is flexibility. By generalizing the sur- 
face definitions and making no assumptions of sym- 
metric mass distributions, the program can be 
applied to any aircraft in an< ase of flight, 
except hover. The unique trimming capababl 1 ity , 
provided by means of a user-supplied subroutine, 
allow unlimited possibilities of trimming strate- 
gies and surface scheduling which are particularly 
important for oblique-winged vehicles and aircraft 
having multiple surfaces affecting a single axis. 
The formulation of the equations of motion permit 
the inclusion of thrust-vectoring effects. The 
ability to select, without proqram modification, 
the state, control, and observation variables for 
the linear models, which when combined with the 
large number of observation quantities available, 
allows any analysis problem to be attacked with 
ease. 
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Main program (low Uteraupplled subroutines 



Fig. 1 Program flou diagram shooing oormunicatitn 
uith ueer-eupplied subroutines. 
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Fig. P. Selection of state variables for linear 
mooel . 


Fig. 3 Inputs and outputs to the user. supplied 
subroutine UCNTRL. 
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Fig. 4 Comparioon of non] inoar and linaar timulationt of a higk-parfoimanea 
aircraft for an alsvator doublet ; laval flight; 0.7 Maoh nunbar; at t 0,000 ft 
altitiidi. 
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